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ABSTRACT 

Recent studies indicate that Type la supernovae (SNe la) consist of two groups — a "prompt" 
component whose rates are proportional to the host galaxy star formation rate, whose members have 
broader lightcurves and are intrinsically more luminous, and a "delayed" component whose members 
take several Gyr to explode, have narrower lightcurves, and are intrinsically fainter. As cosmic star 
formation density increases with redshift, the prompt component should begin to dominate. We 
use a two-component model to predict that the average lightcurve width should increase by 6% 
from z = — 1.5. Using data from various searches we find an 8.1%±2.7% increase in average 
lightcurve width for non-subluminous SNe la from z = 0.03 — 1.12, corresponding to an increase 
in the average intrinsic luminosity of 12%. To test whether there is any bias after supernovae are 
corrected for lightcurve shape we use published data to mimic the effect of population evolution and 
find no significant difference in the measured dark energy equation of state parameter, w. However, 
future measurements of time- variable w will require standardization of SN la magnitudes to 2% up to 
z = 1.7, and it is not yet possible to assess whether lightcurve shape correction works at this level of 
precision. Another concern at z = 1.5 is the expected order of magnitude increase in the number of 
SNe la that cannot be calibrated by current methods. 

Subject headings: cosmology: observations — cosmology: cosmological parameters — supernovae: 
general — surveys 



1. INTRODUCTION 

Type la supernovae (SNe la) are the most important 
standardized candles for cosmology, and have been used 
to discover dark energy and the acceler ating universe 
(|Riess et al.lll998HPerlmutter et al.lll999ft . This was fa- 
cilitated by the realization that supernovae with broader 
lightcurves are intrinsically brig hter, wh i le th ose with 
narrow lightcurves are dimmer ([Phillips! 1 1993T ) . Vari- 
ous schemes exist to correct SN la luminosities based 
on their lightcurve shape (e.g. iRiess. Pre ss. & Ki rshnerl 
19961: iTonrv et al.l 120031: [PTieto, Rest, & Suntzeff 120061: 



Guv et al.ll2007tlJha. Riess. fc Kirshnerll2007t ) -^here we 



use the "stretch" method, in which the time axis of a 
template lightcurve is multiplied by a scale factor s to fit 
the data (jPerlmutter et al.lll997h . 

Some properties of SNe la have been found to correlate 
with environment — brighter supernovae with broader 
lightcurv es (high s) tend t o occur in late-type spiral 
galaxies ([Hamuv et al.| [l995). while dimmer, fast declin- 
ing (low s) supernovae are preferentially located in an 
older stellar population, leading to the conclusion that 
the age of the pro genitor syste m is a key variable affecting 
SN la properties (|Howellll2001[ ). The fact that supernovae 
occur at a much higher rate in late type galaxies, and 
that the SN la rate is pro portional to the core-collapse 
rate (M annucci et al.ll2005t ) is another indication that age 
plays an essential role. Follo wing this previous work, 
IScannapieco fc Bildstenl (|2005l ) model SNe la as consist- 
ing of two populations — a "prompt" component whose 
rate is proportional to the star formation rate of the 
host galaxy, and a second "delayed" component whose 
rate is proportional to the stellar mass of the galaxy. 
iSullivan et alj (|2006l hereafter S06) tie all of these results 
together using data from the Supernova Legacy Survey 



(SNLS), finding that slow declining, brighter SNe la come 
from a young population and have a rate proportional 
to star formation on a 0.5 Gyr timescale, while dimmer, 
faster declining supernovae come from a much older pop- 
ulation with a rate proportional to the mass of the host 
galaxy. 

S06 and lMannucci. Delia Valle. fc Panagial (|2006h pre- 
dict that the SNe la whose rates are proportional to 
star formation will start to dominate the total sample 
of SNe as cosmic star formation increases with redshift. 
Since these SNe are intrinsically brighter, the mean lu- 
minosity of the population should increase with red- 
shift. Here we use the t wo-co mponent SN la model of 
IScannapieco fc Bildstenl (|2005l ) and the stretch distribu- 
tion for each component from SNLS data (S06) to quan- 
tify the expected magnitude of this effect. We then com- 
pare the predicted evolution in lightcurve stretch to SN 
distributions from the SNLS and the Higher-z SN Search 
(|Riess et al.ll2007t l. 

2. PREDICTING EVOLUTION 

IScannapieco fc Bildstenl |2005j parameterize the su- 
pernova rate in a galaxy as 

SNR(t) = AM(t) + BM(t), 

where M(t) is the total stellar mass in the galaxy, M 
is the star formation rate, and A and B are constants. 
These authors use the supernova rates in galaxies of dif- 
ferent morphologies and colors to derive values for A and 
B. S06 use an alternate method — they fit galaxy mod- 
els to SNLS u*g'r'i'z' host galaxy photometry to derive 
masses and star formati on rates. Then, using the c osmic 
star formation history of Hopkin s" fc Beacoml (|2006l ). S06 
predict the rate of SNe from each component versus red- 
shift in their Fig. 10. Note that here we adopt the same 
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Fig. 1. — Stretch distributions of SNe from each component 
of the two component model. Top: Prompt (B-component) SNe 
la. Bottom: Delayed (A-component) SNe la. The distributions 
were derived from those in S06, as described in the text. Best-fit 
Gaussians are shown. The A-component Gaussian is centered at 
s = 0.945, with a = 0.077. The B-component Gaussian is centered 
at 1.071 with a = 0.063. 

definition of M(t) as S06 — it is the mass turned into 
stars and does not include mass loss from supernovae. 

The prompt and delayed SN la components have differ- 
ent stretch distributions (S06). To determine the stretch 
of each SN, here we use the SNe from S06, though we fit 
a new lig htcurve template to the data using the SiFTO 
method (|Conlev et al.l l20d7aH bh . Because stretches are 
always defined relative to the s = 1 template lightcurve, 
stretch values should only compared within a publica- 
tion. However, the stretches deri ved here are a pprox - 
imately 4% larger than those in lAstier et al.l (|2006f) . 
largely due to the use of a narrower s = 1 template. 

All SNe from passive galaxies (i.e. those with no mea- 
surable star formation rate) were assigned to the A com- 
ponent. Star forming galaxies have SNe la from both 
components, so the A distribution from passive galaxies 
was scaled by mass and subtracted from the distribution 
of SNe la from star forming galaxies, leaving the B dis- 
tribution (as in S06). The resulting distributions, and 
Gaussian fits are shown in Figure [TJ Note that to con- 
serve the total number of SNe one should add the SNe 
subtracted from the B distribution back to the A dis- 
tribution - this is unnecessary for our purposes because 
the relative heights of the gaussians are normalized as a 
function of redshift in the next step. 

To estimate the expected stretch evolution with red- 
shift, we take the observed A and B distributions and 
scale them to the predicted relative values with redshift 
from Fig. 10 of S06. Increasing cosmic star formation 
with redshift produces a larger fraction of SNe from the 
prompt component. Stellar mass as a function of redshift 
is determined by integrating the star formation history 
from the earliest times, so the total stellar mass, and the 
number of SNe from the A component, decreases with 
increasing redshift. The net result is that in the A + B 
model the mean stretch increases from 0.98 at z = to 
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Note. — Cols. 2-3: the \ 2 per degree of freedom between the 
data and the predictions of the A+B and modified Mannucci (MM) 
models. Cols. 4-5: KS-test probability that the data is drawn from 
each model. Bins with zero counts were assigned an error of 1.15, 
possibly underestimating the \ 2 - The KS test is also imperfect 
because probabilities were derived for a single Gaussian, not the 
sum of two Gaussians as used here. In both cases the two model 
components were fixed by the A and B numbers in S06. 

1.04 at z = 1.5. 

One caveat is that in the A + B model there is no time 
dependence for the A component. SNe la from 10 Gyr 
old progenitors are just as likely as SNe la from 3 Gyr 
old progenitors. If 10 Gyr-old SNe la are actually more 
rare, the A + B model will overpredict the number of A- 
component SNe at z — 0, as they result from stars formed 
during the high star formation rate in the early universe 
(see discussion in S06). As an alternative to the A + B 
model we test ed the two component SN la delay time dis - 
tribution from lMannucci. Delia Valle. &: Panagial (2006), 
which has an exponential decrease in supernova prob- 
ability from the delayed component with time. The 
drawback of this model is that the probability distri- 
bution is somewhat arbitrary. Also, rather than the 
50-50 split between prompt and d e layed SNe chosen by 
iMannucci, Delia Valle. fc Panagial (|2006f ). here we scale 
each component by the A and B values measured by S06. 
This gives similar results to the A + B model, predicting 
a shift in mean stretch from 0.98 to 1.02 from z = — 1.5. 

3. COMPARISON TO OBSERVATIONS 

In Figure [5] we compare the predicted stretch distribu- 
tions from the A+B model to the observed stretch distri- 
butions in three redshift bi ns from the low redshift data 
used bv lAstier et~aLl (|2006h . the SNLS data in S06, and 
the d ata of the higher-z supernova search (|Riess et al.l 
l2007f l. All lightcurves have been refit here using the same 
method. We also tested the data again st the modified 
IMannucci, Delia Valle. fc Panagial (|2006f ) model, but we 
did not find it to be a better predictor of SN evolution 
with redshift (Table [TJ. 

Each survey has different selection effects - the most 
serious for the current study is Malmquist bias, the ten- 
dency to discover only the brightest members of a group 
near the detection limit of a magnitude-limited survey. 
To minimize the effect, for each of the high redshift 
searches we only consider supernovae from a reduced 
volume so that none of the supernovae used are near 
the magnitude limit. The SNLS regularly discovers SNe 
la out to z > 1, but here we use only the subset with 
z < .75, where Malmquist bias is minim al (Astic r~et alJ 
I2006h . Similarly, we only use lRiess et "all ((20071) SNe with 
z < 1 .5, where the autho rs report their sample is com- 
plete (jStrolger et al.l[2C)04h . Lowering the redshift cuttoff 
to z = 1.2 does not change the average stretch for the 
highest-z sample, but it reduces the sample size from 20 
to 13, and thus decreases the significance of the results. 



Evolution in average SN la properties with redshift 
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Fig. 2. — Actual stretch distributions compared to predictions 
from the A + B model. In each case the A + B model evaluated 
at the median redshift of the distribution is shown. Top: SNe la 
from z = 0-0.1 (median z = 0.026; N=50). Middle: S NLS SNe la 
from z = 0.1-0.75 (median z = 0.55; N=99). Bottom: iRiess et alj 
lt2007h SNe la from z = 0.75-1.5 (median z = 1.12; N=20), with s 
errors < 0.2 (requiring s error < 0.1 reduces the sample to 16 and 
increases the average stretch from 1.06 to 1.07). The vertical line 
gives the mean stretch for each distribution. The top axis converts 
a stretch difference into a magnitude difference using a = 1.5. 

As an additional protection against selection bias, we 
only consider SNe with s > 0.7. SNe la with s < 0.7 
are bot h dim and spectoscop ically peculiar, like SN 
1991bg (|Filippenko et al.lll992t ), and have not yet been 
detected at z > 0.2, probably because of a combin ation of 
Malmquist bias and spe ctroscopic selection bias (jHowelll 
l200lt lHowell"eraT][2005h - as redshift increases, and the 
angular size of the host galaxy decreases, and it becomes 
more and more difficult to spectroscopically identify such 
faint supernovae when blended with their bright, often 
elliptical, hosts. This cut removes 3 SNe la from the low- 
z sample [other low-z SNe are already removed because 
we only consider Hubble-flow SNe la. with z > 0.015, to 
be consistent with lAstier et al.l (|2006l )]. 

In all cases we use only SNe la with at least 4 lightcurve 
points, and at least one detection before 10 rest-frame 
days after maximum light in the B-band, so that stretch 
is accurately measured. 

Figure [2] shows that the average observed stretch in- 
creases with redshift, from s — 0.98 ± 0.02 at a median 
redshift of z = 0.03, to s = 1.02 ± 0.01 at z = 0.56, and 
s = 1.06 ± 0.02 at z — 1.12. Simultaneously the per- 
centage of SNe la with s < 0.9 decreases from 24% to 
15% to 1.4%. The KS test gives a 2% probability that 
the lowest and highest redshift sample are drawn from 
the same distribution. The predicted distributions from 
the A + B model are overplotted. The observed trends 
match the predictions of the empirically-based models — 
with increasing redshift fewer low-stretch SNe la are ob- 




FlG. 3. — Cosmological fits done with the lAstier et"aTI (| 200611 
sample (dashed), and with SNe la drawn from the same sample, 
but using only s < 1 SNe la at z < 0.4 and s > 1 SNe la at z > 
0.4 (solid). Combining with Baryon Acoustic Oscillation results 
(Eiscnstcin et al. 2005) produces the thick lines. The results are 
consistent at the one sigma level. 

served, and the mean SN la stretch increases. We find 
the sa me result when th is analysis is r epeated with the 
SALT (|Guv et al.l l2005h and SALT2 (|Guv et all [2001 
lightcurve fitt ers. These results a re also consistent with 
the findings of lAstier et al.1 (|2006l ). that the low-z sample 
h ad an averag e stret ch 97% that of SNLS SNe. 

lAstier et al.l (|2006l ) estimate distances from SNe la us- 
ing 

Hb = Trig — M + a(s — 1) — /?c, 

where /is is the distance modulus, m* B is the peak B- 
band magnitude, c is a color, and M , a and /3 are param- 
eters fit by minimi zing residuals on the Hubble diagram. 
lAstier et al.1 (2006) found a = 1.52, so a drift in average 
stretch of 0.08 ± 0.026 from z = 0.03 to z = 1.12 results 
in a 12% drift in average intrinsic SN la luminosity over 
this redshift range. 

4. EFFECT ON COSMOLOGICAL STUDIES 

Evolution in the SN population will not necessarily bias 
cosmological studies, since SNe are only used in this way 
after correction for lightcurve shape. However, we can no 
longer assume that any deficiencies in lightcurve width 
correction schemes will average out under the assumption 
the distribution of SNe is similar over all redshifts. If 
there is a systematic residual between low stretch and 
high stretch SNe when they are stretch corrected, this 
could cause a bias in the determination of cosmological 
parameters as the population evolves. 

To test an extreme case of evolution, we fit the equa- 
tion of state parameter for dark ene rgy, w and the mat - 
ter density, Qm using the data from lAstier et all (|2006f ). 
(dashed lines in Fig. [3]) and compared it to a fit using the 
same data, but retaining only s < 1 SNe la at z < 0.4 
and only s > 1 SNe la at z > 0.4 (solid lines). The 
strongly evolving subset gives estimates for w and Om 
consistent with the full set, although the errors are larger 
because there are fewer SNe in the subset. 

One worry with an evolving population is that stretch- 
magnitude or color-magnitude relations may evolve, i.e. 
a or (3 derived at one redshift may not be appropriate at 
another. The values derived here for the strong evolution 
subset and the full set are consistent, but again a strong 
test awaits a larger data set. 

Measuring changes in w with time will require much 
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stricter control of potential evolutionary effects. Possible 
biases depend critically on the exact nature of the evo- 
lution, the experimental design, the cosmology, and how 
the time varia ble component of w is parameterized — e.g . 
w/, w\, or w a (jWeller fc Albrechtl[200llKim et al.ll200l 
iLinderi f2006l ) . However, a good rule of thumb is that 
to keep systematic errors significantly below statistical 
errors for a mission such as SNAP, the corrected magni- 
tudes of SNe la should not drift by more than 0.02 mag 
up to z = 1.7 (|Kim et al.ll2004h . Unfortunately there 
are not enough well measured SNe la in the literature 
to determine whether or not stretch correction works to 
this level of precision. After stretch correction, the rms 
scatter of SNe la around the H ubble line is ~ 0.2 mag in 
the best cases (|Guv et al.ll2007h . Therefore - 100 SNe la 
are required in each of several stretch or redshift bins to 
determine whether biases remain at the 0.02 level after 
correction. Such precise tests will soon be possible by 
combining data from large surveys underway at low and 
high redshift. 

5. DISCUSSION 

We have shown that there is some evolution in the aver- 
age lightcurve width, and thus intrinsic luminosity of SNe 
la from z — to z — 1.5, although significant evolution 
is found only over a large redshift baseline. This evolu- 
tion is consistent with our predictions from the A + B 
model - as star formation increases with redshift, the 
broader-lightcurve SNe la associated with a young stel- 
lar population make up an increasingly larger fraction of 
SNe la. 

Though we have taken steps to minimize the effects 
of Malmquist bias, it is possible that residual effects 
play some role in increasing average stretch with redshift. 
However, the net effect on cosmology studies is the same 
no matter the underlying cause. In cither case, there is 
increased pressure on the light curve shape calibration 
method to correct for the evolution in SN la properties 
w ith redshift. 

iPerlmutter et al.l (| 1999( 1 found that SNe la still give 
evidence for an accelerating universe even if SNe la are 
not corrected for stretch. This was possible because the 
difference in a universe with dark energy and one with 
Qm = 1 is large — 0.25 magnitudes at z — 0.7, whereas 
the population evolution seen here implies that average 
SN la magnitudes should increase by 0.07 mag over the 
same redshift range. However, discriminating between 
dark energy models requires much more precise control 
of SN la magnitudes over a larger redshift baseline. 

Most theoretical studies addressing possible SN 
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la evolution have focused on metallicity. Although 
theorists have proposed various mechanisms that 
could conceivably alter the properties of SNe 
la as the a verage metallicity changes wi th cos- 
mic time (iHoflich. Wheeler, fe T hiclcma ml fl99cl 
Lentz et al.1 l2000t ID ommgucz. Hoflich. & Straniero 
20011: iTimmes. Brown, fe Truranl I2003H there is no 
consensus regarding which effects are important or even 
the sign of these effects. Observational studies have 
found no e vidence that metallicity affects the properties 
of SNe la (lHamuv et al.ll2000t ITvanov. Hamuv. fe Pintol 
[20001 : lOallagher et all 120051 : lEllis et al.ll2007f l. Instead, 
it is more likely that age differences between the two 
populations (and thus almost certainly the mass of 
the secondary star) play a role in the evolu tion of 
the o bserved stretch distribution with redshift (Tipwelli 
[20C31 . 

Prompt SNe la are thought to be brighter because they 
produce more 56 Ni. If the Chandrasekhar-mass model 
describes most SNe la, they must then produce less in- 
termediate mass elements, assuming that the amount of 
unburned material is negligible in normal SNe la. We 
therefore predict that high redshift SNe la will have less 
Ca and Si. This is confirmed b y the most intens ive study 
of high redshift SN la spectra ((Ellis et al.ll2007ft . 

One concern raised by these findings is that patholog - 
ical SNe la s uch as SN 2001 ay dHowell fc Nugentir2 004). 
SN 2 002cx dLi et al.l l200l. SN 2002ic (lHamuv eTaLl 
l20M l. and SNLS-03D3bb ([Howell et al.l l2006fl. which do 
not obey typical lightcurve shape correction schemes, are 
associated with star formation. Since star formation den- 
sity increases by a factor of ten from z = to z = 1.5 
(Hop kins fe Beaco"mll2006t l . at high redshift these patho- 
logical supernovae will be an order of magnitude more 
common. Thus the conventional wisdom that all high 
redshift supernovae will have counterparts at low redshift 
([Branch et al.ll200l( ) only holds if sample sizes are much 
larger than those currently used for cosmolo gy. Only 20 
SNe la have published lightcurves at z > 1 ([Riess et"aH 
1200711 Astier et al.ll2006f L and only ~50 SNe la at z < 0.1 
have sufficient data to be cosmologically useful. Fortu- 
nately, thus far it has been possible to identify these out- 
liers so that they do not affect cosmological analyses, 
but future studies requiring increased precision must be 
vigilant of the effects of an evolving SN la population. 
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